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ABSTRACT 


This package is designed to make all elementary calculations 
involved in evaluating 

N 


,(j) 


(t) = ( l a 1 n 1 k (t)) 

i“l * 


(j) 


» j m 1 » • * » » 


fairly easy where the are usual normalized B-splines. In 

addition interpolation and least squares subroutines are provided 
that one may easily generate "meaningful" spline functions. 



0, Basic Spline Package 


This package is designed to solve some of the elementary problems 
associated with splines (eg. spline interpolation and least squares 
fit) . I have attempted to keep the number of technical subroutines to 
a minimum while giving the user great flexibility in determining the 
"type" of spline desired. This flexibility occurs in the freedom the 
user has in his choice of knot sequence and/or his choice of piecewise 
degree. 

The subroutines fall into three basic categories. The first category 
involves computations with a given spline and/or knot sequence. These 
routines deal with general properties of splines and are not tied down to 
any specific task. Included in this category is 
BSPLVN 
BVALUE 
INTERV , 

The second category involves routines which calculate the coefficients of 
splines which perform a specific task (such as interpolation or least 
squares fit) . In this category are 
EQUATE 
INTERP . 

Finally, in the third category is a banded equation solver BNDSLV which 
has nothing directly to do with splines, but which is needed to solve 
specific linear equations which arise in EQUATE and INTERP. 

The reader interested in more sophisticated applications should 
consult [2], where BSPLVN, BVALUE, and INTERV were obtained and where 


I 


certain program fragments of EQUATE and INTERP were sketched. Programs 

which, for example, solve collocation problems of differential equations 
are also located there. 

The appendix contains a listing of the above subroutines. 
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1. Normalized B-splines 

The normalized B-splines are used almost exclusively for calculating 
and evaluating piecewise polynomials. Here we give a very brief description 
of the normalized B-splines and some of their most important properties. 

For more information the reader should consult [l] p and the references 
therein. 

We set 

(1.1) (s— t ) - 


s-t if t<s 


0 if t>s 


Let k be a positive integer which we will fix throughout this section. The 


biinfinite knot sequence { t . } 


( 1 . 2 ) 


e i 2 ‘i-i 




a - lim t. 
i-+— oo 

b « lim t 

i+oo 1 


^■_oo 


will satisfy the following properties: 


(a or b may be infinite) . 


Although we will always be working on a finite interval with a finite 
number of knots, it is convenient, for the purpose of exposition, to work 
with the above biinfinite sequence. 

The normalized B-spline of order k determined by (t^, • . • , t^ + ^} Is 
denoted by ^ and is defined by 

N i,k (t * " ^i+k^i* ^ t i > * ,,,t i+k^ ^ 8 " t >+ * 


( 1 . 3 ) 



where the k-th divided difference operator [ t ± , . . . is acting on 

the variable s. For example if k*l we obtain 

^ t i+l“ t ^+ " ^ t i" t) + 


(1.4) 


[t^,t^^] (s— 1)^_ = 


fc i + l “ fc i 


Since there might be some confusion as to the value of N (t) at t 

i 1 i^ X 

we arbitrarily declare ^ to be right continuous . 

The following properties of the H ± are very useful: 


(1.5) i) N t k (t) > 0 

ii) supp (N lk ) - [t 4 , C 1+k ] 

00 

m) I N .r k (t) * 1 » te(a,b) . 

i=_co 

Le t us suppose that t fc < t fc+1 and < t^ +l , then the set 


is linearly independent on the interval U k ,t N+1 J . Furthermore, the set 
^ N i,k^i-1 ^ orms a basis for the piecewise polynomials of order k 
(degree < k-1) on the interval [t fc , t^] with knots at {t fc+1> . . . .y and 

continuity conditions as prescribed by those knots. That is, if 


(1.7) 


fc k S * j < 


'J+l 


t j+nH-l - ^+1 



5 


If 1 

then the polynomial pieces must join at in a C K_m " fashion. 

j+1 


if there are no knot multiplicities the N. . ’s span the ^[t, ,t. ,] 

1 1 lc Ic N+x 

piecewise polynomials of order k with knots at (t, ,t„} * 

K+l N 
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2. Evaluating Spline Expressions 

In general, one has an expression of the form 
N 

(2.1) l A. N. (t) = s(t) 

1=1 1 


to evaluate in the interval (t^, t^] (recalling that are 

from some given knot sequence) • Due to the local support property of the 

N i,k' S * 1,5 * Ai) if < t < c j + i then 


j 

(2.2) s(t) - l A.N. . (t) . 

i=j-k+l ’ 


Thus in order to evaluate (2.1) efficiently the first step consists of 
finding the index j such that 


(2.3) < t < t 


The subroutine INTERV accomplishes this. 

In evaluating (2.1) and its derivatives the following formula from [1] 
is useful. 


(2.4) 


(t) 


(k-1) 


Oc-j) 


N 

I 

i-l 


(j) 


N l.k-J (t) 


where 


and 


«> - - v. j . . 


(2.5) 
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Notice that evaluating the derivative of a spline written in terms of 
the normalized B-spline basis reduces to evaluating a lower order spline 
in a similar basis. These handy features are combined in the function 
routine BVALUE. 
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3. Computing the Coefficients of the N 

1 , K. 

In the previous section we indicated how to calculate the value of a 
spline or its derivative if one had the coefficients of the representation,, 
In order to obtain the coefficients one generally has to solve a linear 
system of equations whose matrix entries depend on the values of the 
k or its derivatives. For example, if the interpolating spline is 
desired for the function f, then one needs to solve the system 
N 

(3.1) l A ± N 1>k (x^) • f (Xj ) , j-l,...,N. 

Assuming that the x^’s are ordered then (3.1) will be a banded linear 
system. The matrix will have the form 

(3 * 2) “ij “ N i,k (x j ) * 

Thus it would be useful to have a subroutine which generates simultaneously 
the values of all the nonzero normalized B-splines, The subroutine BSPLVN 
accomplishes this task. Further, the subroutine BNDSLV is used to solve 
banded systems of linear equations. 



4. Spline Interpolation 


This section is a user's guide to the subroutine INTERP. The inputs 
to INTERP (T, A, N,K,TAU > F > IFLAG) are 

T - ( ordered ) knot sequence 

N - number of equations 

K - order of the spline 

TAU - (ordered ) interpolation points 

F - external user supplied function , 

The returns are 

A - N coefficients of the , 

I,k 

I FLAG - I FLAG - 0 or 1* 

This routine solves the linear system of equations 
N 

l A(I)N (TAU(J)) - P(TAU(J)), J - 1 N. 

I«1 

If the linear system is singular then IFLAG is set to 1 otherwise it 
is zero. The user must supply (T(I) : I - l i# ..,N+K} and 
(TAU(I) : I * l,...,N} f both non-decreasing sequences with T(I) < T(I+K) 
and TAU (I) < TAU(I+1), 

It can be shown that the linear system is invertible if and only if 
TAU ( I ) e (T(I) ,T(I+K>). If K is even then a fairly stable interpolation 
procedure is obtained by choosing 
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TAU(I) - T(K) + (1-1)/^ , I - 1.....K/2 

TAU(I) = T(K + (I-K/2)) , I - K/2 + 1 N - K/2, 

TAU(I) - T(N+1) - (N-I)A 2 , I - N - K/2 + 1 N 

where 


A _ T(K+1)-T(K) 

a i m — 


and A, 


T(N-H) - T(N) 
K/2 


INTERP calls BSPLVNj INTERV, and BNDSLV. As the subroutines are 
dimensioned now it Is assumed that K < 10. 



5, Spline Least Squares 

This section is a user’s guide to the subroutine EQUATE. The 
inputs to EQUATE (T,N 8 K,LX,X, G,WEIGHT,A) are 

T - ( ordered ) knot sequence 
N - number of equations 
K - order of the spline 
LX - number of data points 
X - (ordered ) x- co-ordinate 
G - Y-values 

WEIGHT - (positive) vector • 

The return is 


A - N coefficients of the N 

I,K 

The returned coefficients satisfy 


LX / N > 

l WEIGHT(J)[ l A(I)N_ v (x(J)) - G(J) 
J-l \l-l I » K 

(l 


v 

LX 

< l WEIGHT (J) 
' J-l 


(x(J)) - G(J) 


for any vector b - (b,,...,b ) . 

1 n 

It is assumed that (X(J) # J-l, . „ . 8 LX) are ordered an further that 
T(K) < X(l) < X(LX) < T(Nfl). This routine calls BSPLVN and BNDSLV. 
the subroutines are dimensioned now it is assumed that K < 10. 
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6. BVALUE 

The function routine BVALUE will return the value of a spline or any 
one of its derivatives (depending on the calling parameters). The inputs 
for BVALUE(T,A,N,K,X,IDERIV) are 


T - ( ordered ) knot sequence 

A - coefficients of the spline 

N - integer such that X e [T(K) # T(N+1)] 

K - order of the spline 

X - evaluation point 

IDERIV - derivative wanted. 

The return is 


BVALUE (T,A,N,K,X 9 IDERIV) » S (IDERIV) (X) 

N 

where S(-) - J A(I)N_ „(•) . This routine will return a value of zero if 
I»1 

X t [T(K) 9 T(N+1)], This routine calls INTERV, 

This routine is based on formulae (2.2) and (2.5) along with the basic 
identity 


t-t 


(6.1) N 1>k (t) 


t i+k -t i 


» (t) . t l+k~ t 


(t) 
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7. INTERV 

The subroutine INTERV will return the subscript ILEFT which satisfies 
T(ILEFT) < X < T( ILEFT + 1) 

if possible. This routine is used for efficient evaluation of a spline as 
in (2.2). More precisely INTERV (XT „ LXT , X , ILEFT , MFLAG) has inputs 

XT ~ (ordered) vector 

LXT - Index of the right hand end point XT (LXT). 

X - evaluation point. 

The returns are ILEFT and MFLAG where MFLAG = 0 means that 1 < ILEFT < LXT - 1 
and XT (ILEFT) < X < XT (ILEFT + 1). This is the usual case. If MFLAG « -1 
then X < XT (1) and ILEFT - 1. If MFLAG - 1 then XT(LXT) < X and ILEFT * LXT. 
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8. BSPLVN 

The subroutine BSPLVN is used to produce (simultaneously) the values 

of the nonzero normalized B-spline N, , at a given point X. This routine 

1,K 

is used here exclusively to fill out the various coefficient matrices. 

Then inputs to BSP LVN (T , JHIGH , INDEX, X fl ILEFT , VNIKX) are 

T - (ordered) knot sequence 
JHIGH - order of the spline 
INDEX - either 1 or 2 
X - evaluation point 

ILEFT - where T (ILEFT) < X < T( ILEFT + 1), 

The return is 

WIKX(J ) = N I L EFr-JHICH + J,JHIGH (X >- J_1 JH1GH 

if INDEX « 1, This is the only option we have used in EQUATE and INTERP. 

If INDXE - 2 , then 

VNIKX(X) - N ILEFT+I _ JItJ , (» 

where J f « max( JHIGH, J+l) and J is a local variable in the subroutine. This 
option is useful in generating simultaneously the values of some or all of 
the derivatives of the nonzero at X. 

Note that before calling BSPLVN one must find ILEFT so that 
T (ILEFT) < X< T(IEFT-fl) . 
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9, BNDSLV 

The subroutine BNDSLV is designed to solve banded systems of linear 
equations „ The input parameters for BNDSLV(A,B,X,N, IROWLN) are 

A - Banded matrix 

B - right hand side 

N - number of equations 

IROWLN - the row length (band width). 


The return is X a vector which solves 
MX = fi , 


where (setting R = (IROWLN + l)/2) 


M(K,J) 


A(K,R + J - K), | J-Kj < R -1 
0 otherwise 


In short, the diagonal of M is stored in the R-th column of A and the rest of 
the entries are placed accordingly. The routine uses elimination with no 
provision for pivoting. 
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ro 

TO 

90 

0008 


I L 3 

~ L X T — l 


— *- — 


— 

0009 





CO 

70 

?\ 

00 10 

' * 2 0 IF 

(X 

.eg. [ ihi n 

_ „ r . .. . 

GO 

TO 

40 

0011 

21 TF 

( X 

;gf . x t ( it on 


GO 

TO 

1 00 


r****' '('or x‘“ . L" . -xttt'h r ) ;“r- r^irr ag^rr^-n^rnr*"' 


0012 30 IS'FP =1 

001 3 “ 31 IHI - T L 0 - - - - - ~ ■ 77 - z ?* 7 ~ ; 

0014 i to = r hi -• r $- fp v 7 >r r ' y ' ’-. r 

*" 0015 TP “f TIT" 717 . 1 ) ' :^n;~TGr 75 :~" :y J S : "' 77 ' 

0016 IF (X .GE. XT(TLOJ) G3V?0\5O -< ^ 

■^001 7 * TSTT P ~ ~I STF***Z 

0018 . co to 3i ••••';. > . 

0019 ' 3 5 I CD - l ; “C-'T — — -- 

0020 IF (X . 17 . XT ( 1 ) ) . 7 V.G 0 TO 9 q r ? 

0021 ; “ - 7 — GO ^n“ 50 ;-••-••••., * 4 — 

C****_MOW X .O F. XMILOI . FIND V pr *^f ' 80 JMD „,,V >,., r 

0023 41 TLO = I HI 

' 00 24 ' IHI ‘"a- TLO ♦ I STEP " V ~ ; r-— - ~ - - 

0025 IF (IH! *GF . L XT ) '-;G 0 TO 45 

0026 ™ IF IX .LT. XT (IHI)) " v ."GO "TO “50 — “ ~ ” * 

0027 I ST FP = !STEP *2 

“ 00 28 ' ' ' : l "T,n -TT 41 ; 

0029 4 5 IF (X .G 5 . X'MLXrn GO T 0 110 

0030 IHI "=:■ L XT - ~ •' 7 7 ' 7 ”. 7 ~ — . - 

C**** MOW X t f TLO) . L c • X .IT. XT(THI) . HARROW THE INTERVAL 

0031 ' ' 50 MIDDLE = (ILO + IHr )/2 --- — - 

_ 003 ? TF (middle ,E 0 . ILO) GO TO 100 

"C ' M'O r E’« TT "T$ ’ A"5SU'TF0 7H A^"h ID 0LT '"=" TLO ' 1 K^’CVS i: ~ 1HT' =" ICO* I 

0033 !F (X .LT. XT(MIOOLF)) GO TO 53 

'"0034 , ItO'-a MIDDLE • - 

0035 GO to 50 

*0036 ; 53 "" IHI “= "'HITmL'E — - 

0037 GO TO 50 

“ ™" X ? '"- ^ ry ~~GFr~TnTTPT 7 T— ATI D PET U PTr 

90 "FLAG a -1 


0033 




^ ; JFOR TR A N I y G_ IF VF L 2 1 _ I MY F P V 

' 00^2 ' __ II T FT = I L 0 

no ^ F L A 0 = 1 

- mcr- = VX T - r - - 

EMC"' — 


_ - ? * T E = 7 4249 09/77/4 1 

<vr T1J3N 



FORTRAN IV G 


MAIN 


Q4 = 74249 


09/27/40 


ooor 

0002 
0003 
000 4 

0005 

0006 


0007 

0008 
'0009" 

0010 
00 IT 
0012 

0013 

0014 


LEVEL 21 

C ' 

ijT=“TT7 : “Rr“7Tr'Rrr xttofrtv j — - — — 

C ALGOL AT? S VALUE 4 T * X* OF MOFRlV*-rH D?RlV* f IV'* n«= S p LUi - PP n M 

OrMrNSIDN T(lj:,.iCD ■ ' ' * 5 - • 

0 I MfT MS I ON AJ ( 20 ) ,0° ( 20) » 0M( 20) 

n V A LUE - 0 . ‘ ' 

K M I D E Ff = K - 1 0 5 P. I V 

~~~rP“TK'MTDI?R~7t F7 “0 ) “ GO ~t 0“ 99 


^wfinD •SXI4-IN"" (K ,K) SUCH THAT ~TI ) '. LE . X"VtT7 T{ m J” 

(OR I * L E . T (I+1> IF 1 ( 1 ) ,L • . rti + 1 ) » T(N+U). 

“ "KHl "S’ X-l - • • — 

CALL nrCR V J T (K ) , N+1 -K M1, X , I , Me LAG ) 

IF TTLAG) 99 » 20 » 9 

"9 if tx .or.- t( rn 

10 IF ( I .eo. K) 

: f < x . e g . t m ) sfoo; ' 10 ■/ ' * 




09iT0',99 ? r “^-T. 


S V . 


7 V 




0015 

0016 
“0017 
"001 8 
~oorr- 

0020 
002 V“ 
0022 
0023"" 

0024 

0025 

0026 
0027 


002B‘ 

0029 

0030 

0031 
"0032 

0033 

0034" 

0035 

0036 


C **•* 01 FF FT FNC F THE COEFFICIENTS * I D^R I V* f I mf S „ 

2 0 I NK. " = ! -K - r-.^r • - — - • - - 

On ?! J = l rK -\h . » 

“ - - tmkpj jMK-t - j- 7 T 

21 , A J ( J ) = A (IMKPJ) 4 a TTr, 


r p--f idfrt v~mr; o ) — 

22 00 23 J = l, IDER IV 

Knj — 

FKMJ a FLOAT(KMJ) 

THT" =” r~ 

no 23 J J” 1 » KM J 
~TM T's“THITI 




V 

' -T 


’OO'-tCrTT):;"- — - 


23 


IH1KMJ = THI - KM J 

A I ( J J ) = { A J { J J +1 ) AJ ( J J ) ) / (7 ( T H 7 ) . - T ( THMK M J ) ) x? ? K V J 


* * COM PUTT VALUE A T INI Till , ? (I FI J ) np • joERIV-TH "OF 9 IVA T IVF7 

GIV C N ITS RELEVANT 8-SPLINE COEFF. I\* A J(1 A j{ Iv) 

“3fO”iF- ci Dt«rrv-Tor^tr ^go -to ~39 * ------ --i - 

IPl r l+l 

OCT 32 J = 1 ♦ KM I OCR - 

IPJ a I + J 

OP CUT ‘ =~ ? ( T PU ) - X — ; 

IP.IM.J s IPl - j- 

"3 2 CTT^T JT~^ -y- TP-lWTl: • 

I Dt P° l = IOERTV + 1 

n*> 3^ i-r^cupi -*/mi - 


§s 

CA 



742 49 


09 / 27/40 


r> 

FORTRAN 

TV 0 L F VFL 

21 


V 

• f 

0039 


no 

33 JJ=1 

e : 

0040 

0041 

33 


Ajmr 
ILO = I 

• 

004? 

3 9 
r 

S V4 LijT 

T V‘ J ( 1 


004 3 
0044 

99 

FMO 



G , 


VALUF 0 A - c = 

,KMJ 

= ( ?.‘3 i"jj ran v 4 j n .n -op < j j y j / ( 

LC - 1 

) 

. • RETURN 


Mnxoj+DRrjj) > 





FORTRAN IV G LEVEL 21 


’.MAIN . ... oar 


. 0001 
0002 

0003 

0004 

0005 

0006 
0007 
0003 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 

0017 

0018 
0019 
0020' 
0021 
0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 


100 

110 

115 

120 


130 

5 

140 

10 

150 


J M p L r C I T P, F A L * 8 ( A - H , 0 - Z J 
DIMENSION T ( ?o f f f A LM 2 0)" , A (20) 
EXTERNAL F 
PEA 9 (5,100) K,N 
FOP MAT (2T 4) 

J = N + K 


R EAD { 5 , 110) (T(T)iI = l,J) .. .. 

*£*0(5,110) ( TA Uf I ) , I = 1 , N ) 

FORMAT ( 10D8. 2 ) 

WRI 75(6,115) K,N f (T( I), I = l,J>, (TAtj( lit T - 1 

FORMAT ( IX, 1 ,<=« , I 4, * N = • , T 4, 1 3 ( 1 X , D8, 2 ) ) 
CAM I NT F R P ( T y A f N i K • 7 Ati f F f I F L A 0 ) 

I F ( I c L A G • E 0 • 1 ) GO rn 10 

w R ! r E ( 6 , 1 2 0 ) (A { f f, 1 = 1 » N ) ' 

FORMAT (1H 'COEF A' * 10( IX » 010,2 ) ) 

DO 5 Jl = l f 4 ~ ‘ ' . 

DO 5 1=1,5 

x * i -3 ; ‘ 


WRI 75(6,130) X, Y 

FORMAT! 1H V* ACTUAL VALUES »,5X,»Xr~ 
Y-BVALUEC7,AfN*K»X»J) 

WRI T F ( 6 , 140) "j* X, Y / 


soio; 


f _* .P.g ^ a * 9 1 5» 5 Xt 1 X= S_D^10,2 l 5X, ' Y 


wri TP( 6 tl 50) 

FORMAT ( 1H , 'MAT R IX S I NGULA R~* ) 


STOP 

FN3 


F a 74266 


t'N) 





2 , 5 Xf * YV * 


=" ' *010.2) 



INFFRP 


FORTRAN IVG LEVEL 21 


DfilE = 742 


0001 

0002 

0003 

000* 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0.012 
0013 
0 01 * 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 

_0023 

002 * 

0025 

0026 
0027 


13 


1* 

15 


16 

30 


99 


. SUB^niJ.r INF l NT e P P l T » A f N ( K f r aiJ ( F t TFLAO) 

IMPLICIT RFAL*8IA-H # 0-Z > 

01 MENS I ON Ot 20,20) f T (20) « M 2 0 ) , “ A U ( 2 0) , P U ' ' Y (20) 
*< Ml = K - 1 
NP2 M K= V+2-K 

K P< M 1 =< +K — 1 

00 30 I = l | N 

00 13 J= 1 , KPKM1 ‘ ' 

01 T ,J)=0. 

CALL IN-FRV(T(k } f NP2MK,TAU(I ), I L F c ~ , MFL A G ) 

tleft=tlef^+kmi 

1 P ( MFL AG ) 99,15,1* 

IF( I *L 7 .N) GO TO 99 

I LH F r * — 

CALL BSPLVNIT ,K ,1, F ALM I ) , ILF FT, DUMMY) 

L = T L E FT-I ' - 

DO 16 J = 1 , K ' ' 

L-i>i ..rv ■ ‘ , 

OC I » L 1 -DUMMY! J) 

I F ( 0( I , K ) « FO • 0 ) GO F 0 99" ~ 7 ‘ 7 

R(I )-F(T AU( I ) ) , 

CALL 3ND S L V ( 0 f B , A » N , K PK M 1. ) ‘ " r ‘ •/’“V 

IFLAG=0 > 7 

RETURN r r -*.— : 

1 FL AG=1 ’ ,7;. 

RETURN ~ 7 --- — — 

END ■ 



15/05/04 



.FORTRAN IV G LEVFL_ 21 M ft I N * Oft- F = 74266 

v_. 0001 .'^IMPLICIT p E AL * B U-H t O-Z> 

0002 dimension t ( j> o ) , talk 26 ) , a do T “ 

0003 FX rp RNAL F 

0004 READ (5, 100) K ♦ M 

0005 100 FnPMAT( 2 I 4 ) 

0006 J = m + k 

0007 RFAD(5,110) (T(I) t I=l f j) 

0008 P F ft 0 ( 5 » llOI CTAlM 1 1 , I = l , N ) ri - 

0009 110 FORMAT! 1008 * 2 ) 

0010 WRTTF(6*115> K i N * C T ( I) » (TA(J( I ) » I = 1 » N ) 

0°11 U5 FORMft“( ix, • ,K = * , 14, »N=* t T4, 13( 1X,D8.2) > 

0012 CALL INTFRP(T,A t N f K,7AU,F,IFLAG> ' ’ “ 

0013 IF( I c L A G • EQ. 1 ) GO rn 10 

0014 WRIT? (6,120) ( ACT) , I = 1 » N ) ■' 

0015 120 FORMAT ( xh «C06F A ' , 10 { l X , DIO .2 ) ) 

0016 DO 5 J 1 = l f 4 ~ 

0017 _ _• DO 5 1 = 1,5 ^ : v.>;. 

r 0018 ■■■’.“■ j=Jt-i ; ■■■■“ 

ppi9 . ; x-i -3 ■“ ; 

0020 ' Y = F ( X ) ' 7 

. 0021 WRITC(5,130) X,Y , 

0022 130 F OR m a r ( 1 H t • ft C f UA L V A L U E S * » 5 XV * X = • , 010*2*5 X7 * Y = 1 * f) 10 « 2 ) 

, 0023 Y-BVAuietr.fAfN.K.Xfji , v , 

f -' • ' 0024 5 WRITF (6,140) J,X,Y ' — 

___ ^0|| FORMA T (1H * DE R= » , T 5 ,5X, » X« ; « rj 0 1 0. 2 , 5X , « Y= * ,010.2) . 

- 7 . : 0027 10 -WPI TF ( 6 ,150) J'^r‘7/.!' "\ \ 

0028 150 FORMAT (1H , • M AT R t X S I NGljL A R~* I r ~ * 

0029 STOP 

1 0030 " FN3 “ 






